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Abstract 

In a series of recent papers we have used an operatorial technique to describe stock 
markets and, in a different context, love affairs and their time evolutions. The 
strategy proposed so far does not allow any dumping effect. In this short note we 
show how, within the same framework, a strictly non periodic or quasi-periodic effect 
can be introduced in the model by describing in some details a linear Alice-Bob love 
relation with damping. 
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I Introduction and motivations 



In a series of recent papers, [TJ [2j [3j H] , we have introduced and used a number operator 
strategy in the description of some simplified models of stock markets. The original 
reason for adopting this strategy was that all the relevant observables of the markets (i.e. 
the quantity which are used to describe the market) assume discrete values. The same 
framework was recently used for a completely different problem, i.e. for the description of 
love affairs, and in particular for a two-actors (Alice and Bob) relation, and a three actors 
love affair, in which Carla contributes with Alice and Bob to create a love triangle, [51 E]- 
In all these applications the dynamics is generated by a single operator, the hamiltonian 
of the system, which is constructed in a very natural way (see below. This, in our opinion, 
is the most appealing aspect of the procedure) to take into account all the interactions 
between the different agents of the model. However, the systems considered in [51 [6] are 
closed: the actors have no connection with the environment! In particular, this means 
that Alice can only interact with Bob (and with Carla), but she cannot interact with 
anyone else. Of course, to be more realistic, this is a point to be reconsidered: Alice's 
LoA (level of attraction, see below) for Bob does not only depend on Bob himself, but 
also on all the other possible interactions that Alice may experience with different people, 
as well as on different states of mind or other psychological effects. The same holds true 
for Bob, obviously. In this short note we show how these extra effects can be taken into 
account within the same general framework we have used so far, and how the presence 
of the environment introduce the possibility of getting solutions which go down or up in 
time to some value fixed by the environment itself. In order to make the details of our 
approach completely explicit we will discuss here a linear model for which the solution 
can be deduced analytically, while we postpone the analysis of a non linear model to a 
paper which is now in preparation, where we will use numerical computations, [7]. 

Of course one may wonder why we should use quantum mechanics in modeling classical 
systems. The answer, at least as far as our approach is concerned, is that it is quite easy 
to write down a single operator, the hamiltonian of the system, which fixes the dynamics 
of the system: it is enough to follow few simple rules, [1]- [6]. We will come back to this 
aspect later on. We also want to stress that in recent years quantum mechanical tools 
have been used more and more in connection with many different classical (complex) 
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systems. We cite here just two recent monographs, [5J H], where other references can be 
found. More references can be found in 

The paper is organized as follows: 

In Section II, after a short review of some results of [5], useful to fix the notation, we 
introduce a new model in which the two actors, Alice and Bob, interact among themselves 
and with their own environments. We deduce some of the features of this model, in 
particular the existence of a constant of motion, and we obtain the related equations of 
motion, working in the Heisenberg picture. We also find the explicit solution for these 
equations. 

In Section III we comment these results. Then we give our conclusions and we discuss 
our plans for the future. 

To keep the paper self-contained, we give in the Appendix some basic facts on quan- 
tum mechanics and we discuss how an exponential law can be derived by rather general 
assumptions. Of course, this Appendix is devoted only to those readers which are not 
familiar with these arguments. 

II Old and new models 

In [5] we have introduced a simple model of two lovers, Bob and Alice, who interact 
exhibiting a certain interest for each other. Of course, there are several degrees of possible 
interest, and to a given Bob's interest for Alice (LoA, level of attraction) , there corresponds 
a related reaction (i.e., a second LoA) of Alice for Bob. 

Using a slightly different, and more convenient, notation with respect to that adopted 
in [5], we now introduce a and b, two independent bosonic operators. This means that 
they obey the commutation rules 

[a,<f] = [b,b J t] = l, (2.1) 

where 1 is the identity operator, while all the other commutators are trivial. Recall that, 
for two operators x and y, [x, y] = xy — yx. Here a stands for (the annihilation operator 
associated to) Alice and b for (the annihilation operator associated to) Bob. As in [TJ- [6] , 
bosonic operators are useful since, see Appendix, they create or annihilate excitations or, 
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in our language, modify the LoA's of both Alice and Bob. Further, let y?o be the vacuum 
of a and b, aipo = b(p = 0. If the system is described by (po, then Bob doesn't experience 
any attraction for Alice and viceversa, see below. Using ipo and the creation operators a) 
and b' we may construct the following vectors: 

^^-^(bW^, (2.2) 

where rij = 0, 1, 2, . . ., and j = a, b. Let us also define N a = a, N b = b, and N = 
N a + N b . Hence (see the Appendix), for j = a, b, NjLp najnb = nj(p na>rib , Nj Oj (p na ,n b = {rij - 
1) a jVn a ,n b , and Nj a] (p na ,n b = (rij + 1) a]cp na>nb . We also have N<p nainb = (n a + n b )cp najnb . 
As usual, the Hilbert space % on which the operators act is obtained by taking the closure 
of the linear span of all these vectors, for rij > 0, j = a, b. A state over the system is a 
normalized linear functional oo na .n b labeled by two quantum numbers n a and n b such that 
UnaiUbiz) = (fn a ,n b ,xip najnb ), where (., .) is the scalar product in 7i and x is an arbitrary 
operator on H. Of course, for generic x, u n , a ^ nb (x) is a complex number; if x — x* , which 
is what happens in all our computations, u na! n b (x) is a real quantity. 

In this paper we associate the (integer) eigenvalue n a of N a to the LoA that Alice 
experiences for Bob: the higher the value of n a the more Alice desires Bob. For instance, 
if n a = 0, Alice just does not care about Bob. On the other hand, we use n b , the eigenvalue 
of N b , to label the attraction of Bob for Alice. In [5] we have assumed that the dynamics 
of the love affair can be deduced by the following hamiltonian 

H = A (a M tf + h.c.) , (2.3) 

where A is the interaction parameter, which also play the role of a time scaling, [5]. 
Notice that no free hamiltonian is considered here: H consists only of the interaction 
contribution. The physical meaning of H can be deduced considering the action of, say, 
a M b^ on the vector describing the system at time t — 0, (p na ,n b - This means that, at 
t — 0, Bob is in the state n b , i.e. n b is Bob's LoA, while Alice is in the state n a . However, 
because of the definition of (p na ,n b i aM ^ (p na ,n b , which is different from zero only if M < n a , 
is proportional to (p na -M,n b +i- Hence, Alice's interest for Bob decreases of M units while 
Bob's interest for Alice increases of 1 unit. Of course, because of the presence of ba) M in 
H, if n b > 1 we see that ba) M ip na ,n b is proportional to ip na+ M,n b -i' hence, Alice's interest 
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for Bob is increasing (of M units) while Bob looses interest in Alice. Hence H describes 
quite well what in [S] has been assumed to be a natural law in love affairs: the more Bob is 
attracted by Alice, the less Alice cares about Bob, and viceversa. I(t) := N a (t) + M Ni(t) 
is a constant of motion: I(t) = 1(0) = N a (0) + M N b (0), for all t > 0: [H, I] = 0. 
Therefore, during the time evolution, a certain global attraction is preserved and can only 
be exchanged between Alice and Bob. Notice that the existence of such a constant of 
motion is in agreement with our naive vision of a love affair sketched above. 

The hamiltonian (12. 3p produces the following differential system for the annihilation 
operators a(t) and b(t): 

a(t) = -i\Mb(t)(a(ty) M - 1 f24l 
bit) = -i\(a(t)) M . 1 ' ' 

Then, we may use the solutions of these equations to construct N a (t) := a^(t)a(t) and 
Nt,(t) := b^(t) b(t). The operators a(t) and b(t) can be found analytically if M — 1. In this 
case, which corresponds to the assumption that Alice and Bob react in the same identical 
way (but for a sign), the solution of system (12 .4p is 

a(t) = acos(Xt) — ibsin(Xt), b(t) = bcos(Xt) — iasm(Xt). (2.5) 

Now, if we assume that at t = Bob and Alice are respectively in the n a 'th and n&'th 
LoA's, the state of the system at t = is oj najnb . Therefore, calling rij{t) := u nat n b (Nj{t)), 
j = a, b, we find that 

n a (t) = n a cos 2 (\t) + rib sin 2 (\t), n\,(t) = rib cos 2 (At) + n a sin 2 (At). (2.6) 

Hence u na:nb (I(t)) = n a + n b , as expected. The conclusion is quite simple and close to 
our view of how the law of the attraction should work: the infatuations of Alice and 
Bob oscillate in such a way that when Bob's LoA increases, that of Alice decreases and 
viceversa, with a period which is directly related to the value of the interaction parameter 
A, which therefore can be seen as a time scaling. In particular, as it is natural, if A = 
equation (12.61) shows that both Alice and Bob stay in their initial LoA's. 

Much harder is the situation when M > 1, [5], for which a numerical approach seems 
more appropriate. Here we are interested in introducing within our scheme the role 
of the environment of both Alice and Bob, in a such a way that the two main actors 
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of our game can get in touch with the external world. We believe that, doing so, the 
model approaches more and more a real system. However, since here we are interested in 
describing the general ideas of our method, we will only consider linear interactions, i.e. 
we fix M = 1. In this way an analytical solution will be found and the general strategy 
will appear clearly with no extra difficulties arising from non linearity. In a paper in 
preparation, [7], non linear effect will also be considered and a numerical approach will 
be proposed. 

II. 1 Enriching the model 

Suppose now that Alice and Bob interact with their own reservoirs. We call S = S a U 
Sb U IZa U7^b the full system, made of Alice (S a ), Bob (<S&), Alice's reservoir (Ha), and 
Bob's reservoir (TZb)- The hamiltonian for S looks now like 

H = Ha + Hb + A-ffj, 

H A = Lo a a^a + J R Vt A (k)A\k)A(k) dk + ^ A J R (a)A(k) + aA\k)) dk, 
H B = uJbtfb + f M n B (k)B\k)B(k) dk + lB f M (b^B(k) + bB\k)) dk, 
Hi = a)b + atf. 

All the constant in (12.71) are real quantities. It is worth mentioning that a regularization 
could be also considered in the definition above to make the hamiltonian rigorously de- 
fined. We will skip this mathematical details here, since we are more interested in the 
physical meaning of H and since, using for instance the stochastic limit approach, [TU], we 
can make our treatment rigorous. The following bosonic commutation rules are assumed: 

[a, at] = [b, 6+] = 1, [A(k), A\q)} = [B(k), B\q)) = 15(k - q), (2.8) 

while all the other commutators are zero. Ha and Hb respectively describe the interaction 
of Alice and Bob with their own reservoirs, which consist of several (infinite) ingredients. 
Their forms are, in a sense, standard for systems interacting with a reservoir, see [12] for 
instance. These can be interpreted as other men, other women, or as the effect of the 
tiredness or tediousness in the sentimental relation, as well as some other (psychological) 
effects. Hi contains our (linear) interaction between Alice and Bob, which follows the 
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same rules of our previous model: for instance, a) b describes an increasing Bob's LoA 
and a simultaneously decreasing Alice's LoA. 

The Heisenberg equations of motion for the annihilation operators are the following 



a(t) = i[H, a(t)] = —iu a a(t) — i^ A J R A(k, t) dk — i\b(t), 
b{t) = i[H, b(t)} = -iu b b(t) - J R B(k, t) dk - i\a(t), 
A(k, t) = i[H, A(k, t)] = -iQ A (k)A(k, t) - i-f A a(t), 
[ B(k, t) = i[H, B(k, t)} = -iQ B {k)B(k, t) - i^ B b{t). 



(2.9) 



Using the same strategy discussed in Appendix, i.e. choosing Q A (k) = &>Ak and Q B (k) 
Q B k, fl A , > 0, we further get 




-v A a(t) - i\b(t) - i^ A f A {t), 
-u B b(t) - i\a(t) - ij B f B {t), 



(2.10) 



where v A = v B = f A (t) = J R e-^ kt A(k) dk and f B (t) = ^e~^ kt B{k) dk. 
From (I2.10p we deduce the following second order differential equation for a(t): 

a(t) + d{t){u A + v B ) + a{t){u A u A + A 2 ) = $(t), (2.11) 

where $(t) := — i'jAfAf) — ^bIaJa^) ~ -^7b/s(*)- To find the solution of this equation, 
satisfying a(0) = a and 6(0) = b, we need to compute the Green function for the differential 
operators L := -C- + (v A + v B ) 4 + [y A v A + A 2 ), i.e. the function G(t) satisfying L[G](t) = 
S(t). To keep the computations reasonably simple from now on we fix v A = v B =: v, 
which implies that = and oo A = lo b =: u>. It is a standard exercise in Fourier 
transform to deduce that 



G(t) 



i sin(At) e~ v \ t > 

0, otherwise. 



The general solution of the equation L[a (t)] = is a>o(t) = x + e e+t + x_e e ~*, with e± = 
{—v ± iX). After few computations and recalling that b(t) can be deduced from a(t) by 
rewriting (I2.10p as b(t) = j- (ia(t) + iv A a{t) — j A f A (t)), we get 

a(t) = ae~ ut cos(At) - ibe~ ut sin(At) + R a (t), 
b(t) = be~ ut cos(At) - iae~ ut sin(At) + R b (t), 
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where we have defined the following functions: 



' R a (t) = pit) + 2=^ {zr(0) sin(At) - Ap(0) cos(At)} , 
R b (t) = ± {T(t) - T(0)e- ut cos(At)} + tpiO)e~ ut sin(At), 
pit) = J R (p A (k,t)A(k) + p B (k,t)B(k)) dk, 

k r(*) = / R (r A (fc,t)A(fc) + r fl (fc,t)5(fc)) dfc, 



(2.13) 



as well as 




p B (k,t) 




[ T A (k,t) =ip A (k,t) +iup A (k,t) --f A e- iQAkt , T B (k,t) = ip B (k,t) + iup B (k,t). 

It is now easy to find the mean value of the number operators N a (t) = a'(t)a(t) and 
Nf)(t) = tf(t)b(t) on a state over S, which is of the form cus(-) '■= (<Pn a ,n b , ■ Vn a ,n b ) ^ni-)- 
Here ujk is a suitable state over the reservoir 1Z = 1Z A U TZ B , while, as usual, <f na ,n b is the 
eigenstate of N a and iV& with eigenvalues n a and n&. Then, calling n a (t) = us(N a it)) and 
rib{t) = us(Nb(t)), and assuming that, see Appendix, 



u n {A\k)Aiq)) = N A ik)Sik - q), u n iB\k)Biq)) = N B (k)6(k - q), (2.14) 



N A ik) and N B ik) to be fixed, we conclude that 



n a (t) = e-^/nx (n aCOS 2 




2 (At)) + 




(2.15) 



and 




sin 2 (At) + 



(2.16) 



where we have introduced the following functions: 



0o,a(M) = MM) + ^e- l/t sin(At)r j4 (fc,0) 

Mo,b(A:, t) = p B (A:, t) + % -e~ vt sin(At)r B (A;, 0) 

A 

fi>b,A{k, t) = j T A (k, t) - je~ vt cos(At)r A (A;, 0) 
W,b(M) = T T B (k,t) - \e- ut cos(Xt)T B {k,0) 



e~ vt cos(At)p B (A;,0) 



e~ vt cos(At)p A (fc,0), 



e _iyt sin(At)/u(fc,0) 



e - ^ sin(At)p B (A;,0) 
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Formulas (I2.15P and (I2.16P look interesting since they clearly display the different contri- 
butions arising from the system S a U S b and from the reservoir 1Z. Of course, different 
choices of the functions of the reservoir NA(k) and Ns{k) clearly produces different ex- 
pressions for the LoA's of Alice and Bob. In particular, if for instance A^/c) = Ns{k) = 
almost everywhere in A;, we get 

n a (t) = e~ 27T ^ t/nA (n a cos 2 (\t)+n b sm 2 {Xt)) , 

n b (t) = e~ 2 ^ t/QA (n b cos 2 (At) + n a sin 2 (At)) , (2.17) 

which produce damped oscillations for both Alice and Bob: independently of their initial 
status, the effect of the reservoirs is to switch off the love between Alice and Bob, at least 
for this trivial choice of Aa(/c) and Ns{k). In this case, if for instance, n a = n b , i.e. if 
Alice and Bob experience the same LoA at t — 0, then we get n a (t) = n a e -271 " 7 ^/^ anc [ 
n b (t) = n b e~ 27r TA*/^-4 _ ^he S p eec l Q f decay of their LoA is related to j a /Qa which, we 
recall, coincides with 7^ /fig. In particular, the stronger the interaction between, say, 
Alice and her reservoir, the faster the decay to zero of her love for Bob. Of course, a 
different speed (i.e. behavior) for Alice and Bob is expected, in general. This can be 
recovered assuming that ^ ^ vb- This situation will be considered in [7]. 

Ill Final comments 

Let us now consider what happens if N^{k) and Ns{k) are different from zero. First we 
consider the following case: A^(fc) = Ae(fc) = 5e~ fc2 . Here 5 is chosen because, just to 
fix the ideas, we will consider in the following that, at t = 0, the initial conditions on n a 
and n b are simply 1 or 5: n a — 1 and n b = 5, see Figure (U corresponds to a low Alice's 
LoA and an high Bob's LoA. The other parameters are chosen to be the following: u — 1, 
Qa — = 1) Ja — 7b = -1 an d A = .3. So the situation looks rather symmetrical and 
the only asymmetry is given by the initial conditions above. 

These figures show a non purely oscillatory behaviors of n a (t) and n b (t), which seem to 
converge to a limiting asymptotic value for large t: this appears to be a sort of intermediate 
state (a sort of attractor) in which both Bob and Alice will continue their relationship 
with not many changes in their LoA'a. However, different situations may arise. First of 
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10 20 30 40 50 10 20 30 40 50 

Figure 1: n a (t) (left) and n b (t) (right), for n a = 1 and n b = 5 



all, if we just change the initial conditions requiring that n a = rib — 5, it is easy to see 
that n a (t) = n b (t) and that they both decay monotonically to a value close to 2.5. 

Moreover, if we now break down the original symmetry between Alice and Bob not only 
considering different initial conditions but also changing the values of the parameters, then 
the asymptotic behavior of n a (t) and rib{t) is not so clear, see Figure [2J This is obtained 
taking u and A as above, and Qa — 1, — 4, ja — -1 an d 7b = .2. This choice satisfies 
the requirement ^\/VLa = IbI^b- 




10 20 30 40 50 10 20 30 40 50 

Figure 2: n a (t) (left) and n b (t) (right), for n a = 5 and n b = 1 



We see from the figure that n a (t) and n b (t) do not apparently tends to a given value for 
t increasing. Similar plots are obtained if we increase even more the asymmetry between 
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Alice and Bob by considering their interactions with two different reservoirs. Indeed, if we 
take now A^i(fc) = 5e _fc2 and Nsik) = 5e~ 3fc2 , leaving all the other constants unchanged, 
we get the functions plotted in Figure EJ for n a — l and = 5, in Figure IU for n a = 5 
and rib — 1) an d in Figure El for n a = 5 and rib — 5. 

In particular Figure [5] shows that n a (t) and n^t) both decrease, not monotonically, 
and that some limiting point, if any, is reached very far away (compared with Figured]) in 
time. Moreover, this same figure also clearly displays that n a (t) +n&(t) is not a constant 
of motion, as expected since the contribution of the reservoirs cannot be neglected. In 
fact, we can show that an integral of motion still exists for this open system, and it looks 
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10 20 30 40 50 10 20 30 40 50 

Figure 5: n a (t) (left) and n b (t) (right), for n a = 5 and n b = 5 

like J = N a + N b + N A + N B , where N a = at a, N b = &t 6, jV^ = J R At(fc)A(fc) dfc and 
Nb = J R B'(k)B(k) dk: [J, H] = 0. This implies that the oscillations and the decay of the 
functions n a (t) and n b (t) must be related to some oscillations and to an overall increasing 
behavior of N^t) and Ns(t). Therefore, Alice's LoA decreases because she interacts with 
other sources. The same happens to Bob. However, if they are lucky enough, their LoA's 
both oscillate around some strictly positive value, as in Figure [T] for instance. So their 
love relation can survive for a long time, even if with less strength that at the beginning. 
This is realistic, indeed. If they are not lucky, see formulas ( I2.17p . their love just vanishes 
after some oscillations. Our analysis shows that, to be concrete, their luck really depends 
on the parameters of their reservoirs. 

In this paper we have continued our analysis of sentimental relationships using oper- 
atorial techniques, first discussed in [5]. In our opinion, our previous models have been 
made significantly more realistic by adding the possibility of damping effects which were 
not there in the original scheme. These effects are related, as we have seen, to the possi- 
bilities for Alice and Bob to interact with the outer world, i.e. with friends, other lovers, 
parents and relatives for instance, which in a realistic love affair play a relevant role in 
the evolution of the relation. 

This improvement will be used in the forthcoming construction of even more realis- 
tic (probably non-linear) models and in the analysis of other biological and economical 
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systems, for which, however, producing exactly solvable models is much more difficult. 
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Appendix: Some general results 
A.l Mathematical background 

Let % be an Hilbert space, S our physical system and 21 the set of all the operators 
on % useful for a complete description of S (the observables of S). In many physical 
systems 21 consists also of unbounded operators. This creates mathematical difficulties 
which will not be considered here, since the relevant part of their spectra are bounded. 
The description of the time evolution of S is related to a self-adjoint operator H = W, 
which will be assumed not to depend explicitly on time, which is called the hamiltonian 
of S. In the Heisenberg representation the time evolution of an observable X £ 21 is 
given by X(t) = e lHt Xe~ lHt or, equivalently, by the solution of the differential equation 
MM = ie iHt [H,X}e- iHt = i[H,X(t)], where [A, B] := AB - BA is the commutator 
between A and B. 

In this paper a special role is played by the so called canonical commutation relations 
(CCR): we say that a set of operators {ai, a], I = 1,2, ... ,L} satisfy the CCR if the 
following hold: 

[a h ajj = S in l, [a h a n ] = [aj, ajj = 0, (A.l) 

for all l,n = 1,2, ... ,L, see [11]. The operators hi = ajai and iV = ^2f =1 ni are both 
self-adjoint operators. In particular hi is the number operator for the 1-th mode, while iV 
is the number operator of S. 
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The Hilbert space of our system is constructed as follows: we introduce the vacuum 
of the theory, that is a vector tp which is annihiled by all the annihilation operators af. 
aitfo = 0, for I — 1, 2, . . . , L. Then we act on <po with the creation operators a\: 

<p ni , n2 ,..., nL := ) . {a\) ni {a\r ■ ■ ■ (4)"^o (A.2) 

These vectors form an orthonormal set and are eigenstates of both hi and N: hiip ni>n2t ___ inL = 
nip ni ,n 2 ,...,n L and Nipn U n 2 ,...,n L = Ntp ni >n2 ,...,n L , where iV = Ya=i u i- Moreover we have 
^(°Wni,n 2 ,...,nJ = (ni ~ l)(a;W,n 2 ,...,nJ and rzj(<Wni,n 2 ,...,nJ = {n L + l){a\y nun2 ,,.,, nij ). 
The Hilbert space is obtained by taking the closure of the linear span of all these vectors. 

An operator Z e 21 is a constant of motion if it commutes with H . Indeed in this case 
the equation of motion for Z(t) reduces to Z(t) = 0, so that Zit) = Z for all t. 

The vector tp nijn2) ... jnL in (|A.2[) defines a vector (or number) state over the algebra 21 

as 

where < , > is the scalar product in We refer to pQ for further details on this subject. 



A.2 An exponential law 

Suppose now that we have a first system, S, interacting with a second system, S, and 
suppose that both S and S are of the same size: by this we mean that S describes a 
single actor whose related (bosonic) operators are a, and h a = a^a and, analogously, 
S describes a second actor whose related (again, bosonic) operators are b, and hb = 
b^b. These operators obey the following rules: [a, a)] = [b, &' ] = 1, while all the other 
commutators are zero. A natural choice for the hamiltonian of the two coupled systems 
is h — u a h a + Ubh b + fi {a)b + b^a) , where u a , u b and /i are real quantities, in order to have 
h = h' . h contains a free part plus an interaction which is such that, if the eigenvalue of h a 
increases of one unit during the time evolution, then the eigenvalue of n& must decreases 
of one unit, and viceversa. Moreover [h,n a + h^] = 0, so that h a + is an integral of 
motion. The equations of motion for a(t) and b(t) can be deduced as follows 

a(t) = i[h,a(t)] = —iu a a(t) — i/ib(t), b(t) = i[h,b(t)] = —iuj^bit) — ifia{t), 
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whose solution can be written as a(t) = a a (t) a + a b (t)b and b(t) = (3 a (t) a + (3 b (t)b, 
where the functions <x,(t) and /3j(t), j = a,b, are linear combinations of e A± *, with A± = 
f{u a + cu b - y/(cu a - cu b ) 2 + V). Moreover a o (0) = (3 b (0) = 1 and a b (0) = (3 a (0) = 0, so 
that a(0) = a and 6(0) = 0. Hence we see that both a(t) and b(t), and n a (t) and n b (t) as 
a consequence, are linear combinations of oscillating functions,: no damping is possible 
within this simple model. 

Suppose now that the system S is replaced by an infinitely extended reservoir TZ, whose 
actors are described by an infinite set of bosonic operators b(k), b^{k) and n(k) = b^(k)b(k), 
fc6M. The hamiltonian of S plus 1Z which extends h above is 

H = H + \H I , H = uh a + I cu(k)n(k)dk, Ej = [ (atf{k) + a*b(k))f(k)dk, 

(A.4) 

where [a, a/] = 1, [b(k) , (q)] = \5{k — q), while all the other commutators are zero. All 
the constant appearing in ( 1A.4[) . as well as the regularizing function f(k), are real, so that 
H = '. At a certain point in this appendix we will take f(k) = 1: this will make easier 
the computation of the solution of the equations of motion, but makes our H a formal 
object. A more rigorous approach, which we will not consider here, can be settled up using 
the results in [10]. Notice that an integral of motion again exists, 9T := n a + J R n(k) dk, 
which extends the previous one. 

With this choice of H the equation of motions are 

d(t) = i[H, a(t)} = —iua(t) — i\ J R f{k) b(k, t) dk, 
b(k,t) =i[H,b(k,t)] = -iu{k)b{k,t) -iXf(k)a(t), 

which are supplemented by the initial conditions a(0) = a and b(k, 0) = b(k). The second 
equation in ([Oil can be rewritten as b(k,t) = b(k)e~ iujik)t -i\f(k) Jja^ije-^^^J dt x . 
Fixing now f(k) = 1, assuming that cu(k) = k, and replacing b(k,t) in the first equation 
m flA~5l) . we find that 

a{t) = -{iu + n\ 2 )a{t) -iX I b{k) e~ ikt dk. 

Jr 

Here we have also changed the order of integration and we have used the following equal- 
ities: J R e~' lk ^~ tl ^ dk = 2nS(t — t\) and Cg(t\)8{t — ti) dt\ = \g(t), for any test function 
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g(t). Then, its solution can be written as 

a(t) = (a-i\ [ dkr){k,t)b(k) ) e ~(^+^' 2 )\ (A.6) 



where r)(k,t) = ^e^ 1 — l) and p(k) = i(u — k) + nX 2 . Using complex contour 
integration it is possible to check that [a(t), a^(t)] = 1 for all t: this means that the natural 
decay of a(t) described in ( 1A.6I) is balanced by an opposite reservoir contribution. This is 



expected because of the existence of the integral of motion 91, and it is crucial since it is 
a measure of the fact that the time evolution is unitarily implemented as it must be since 
H is self-adjoint, even if a(t) goes to zero with t. Let us now consider a state over S (g) 7Z, 
(X s <g) X-ji) = (ip na , X s ip na ) (Xn) n , in which X s and X-u are, respectively, operators of 
the system and of the reservoir, (f Ha is the number eigenstate of n a , and < >^ is a state of 
the reservoir, which is assumed to satisfy ( b\k)b{q]) n = rib(k)5(k — q). This is a standard 
choice, see for instance [12], which naturally extends to TZ the choice we have made for S. 
Then, if we take n b {k) to be constant in k we get, calling n a (t) =< h a (t) >=< a^(t)a(t) >, 

2X 2 Trt i „ ( 1 -2X 2 nt 



n a (t) = n a e- 2 ™ + n b [I - e" 2A ^ j , (A.7) 

which goes to n b as t — > oo. Hence, if < n b < n a , the value of n a (t) decreases with 
time. If, on the other way, n b > n a , then the value of n a (t) increases for large t. This is 
the exponential rule which, as we have shown before, cannot be deduced if TZ has not an 
infinite number of degrees of freedom. 

It might be interesting to remark that the continuous reservoir could be replaced 
by a discrete one, in which we have again an infinite number of actors, but they are 
labeled by a discrete index. In this case, to obtain a Dirac delta, rather than the integral 

use the Poisson summation formula, which we write 
here as ^2 n€ z,e mxc = Y2 n &z^ { x ~ n lf)- However, we will not consider this possibility 
here. 
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